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ABSTRACT 

An iterative method for blade design based on an Euler solver and described in an earlier 
paper is used to design compressor and turbine blades providing shock free transonic flows. 
The method shows a rapid convergence, and indicates how much the flow 7 is sensitive to small 
modifications of the blade geometry, that the classical iterative use of analysis methods might 
not be able to define. 

The relationship between the required Mach number distribution and the resulting geometry 
is discussed. Examples show 7 how 7 geometrical constraints imposed upon the blade shape can 
be respected by using free geometrical parameters or by relaxing the required Mach number 
distribution. 

The same code is used both for the design of the required geometry and for the off-design 
calculations. Examples illustrate the difficulty of designing blade shapes with optimal perfor- 
mance also outside of the design point, 
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INTRODUCTION 

The design of new compressor and turbine blades is in most cases still done by successive 
direct analysis of the flow field around a given blade shape and modifications of the blade 
geometry, according to some empirical criteria and/or the designer’s own experience. This 
approach makes it easier to respect geometrical and mechanical constraints imposed to the 
designer, such as thickness distribution, inertia momentum, stagger angle, pitch-to-chord ratio, 
etc. 

New aerodynamic design tools have been developed, that have shown the ability to provide 
conclusive improvements of the aerodynamic performance when compared to existing blades. 
These improvements result from a specified controlled diffusion along the blade surface or 
a shock-free transonic flow. It is unlikely that they can be obtained by a traditional design 
procedure, namely by a series of flow analysis and empirical blade modifications. The design of 
transonic shock-free blades by means of an inverse method is one of the main topics discussed 
in this paper. 

Analytical design methods developed in the past, using conformal mapping (Lighthill, 1945, 
Woods, 1955), permitted to build a complete theory of the inverse design of airfoils and blades, 
and provided the conditions required for the existence of a solution. However, they have a 
limited application due to the restrictive assumptions needed to allow an analytical formulation 
of the problem. As the blade shape results from the calculation, it is also more difficult to 
satisfy the mechanical constraints that one may wish to impose on the blade shape. 

Numerical inverse methods have been developed for potential flows, using singularities for 
incompressible cases (Murugesan and Railly, 1969, Ubaldi, 1984, Van den Braembussche et 
al., 1989) and the odograph plane (Bauer et al., 1972, Sanz, 1984) or the potential-stream 
function plane (Stanitz, 1953, Schmidt, 1980) for the compressible cases. The last methods 
are not very accurate in the stagnation point region and are unable to predict shocks. It is 
therefore questionable whether blades designed in this way for shock-free transonic flows are 
shock-free in reality. 

Non potential flow fields require solving the Euler equations. Such methods are capable of 
treating shocks correctly and are therefore suited to verify shock-free designs. They are mostly 
used in iterative procedures and require a first guess of the blade shape. This initial geometry 
is modified from the results of a flow analysis until the imposed pressure or velocity distribution 
is reached. The blade modifications can be calculated in a pure mathematical way, in order to 
minimize an error function, eg. depending on the difference between the calculated pressure 
distribution and the target (Vanderplaats, 1979, Hicks, 1981). Although these methods have 
the capability to respect geometrical constraints, they are still very expensive in terms of CPU 
time, because many iterations and flow analyses are required. 

The blade modification can be determined in a more physical way, resulting in decrease in 
CPU time. The present method imposes the required Mach number distribution as a boundary 
condition on the blade wall and uses the concept of a permeable wall tc define the modification 
of the geometry. This approach allows a reduction of the number of blade modifications, and 
consequently of the number of mesh generations. The method has proven to be very efficient 
in subsonic and transonic applications (Leonard, 1990, Leonard and Van den Braembussche, 
1991). As shown in this paper, the iterative procedure makes it easier to meet geometrical and 
^mechanical constraints imposed in industrial applications, and to find out whether a realistic 
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blade shape corresponds to the required Mach number distribution. Another advantage of the | 
present method is the possibility of using the same code for the blade design procedure as well 
as for the off-design analysis. 

THE EULER SOLVER 

The system of Euler equations for unsteady flows is solved using a time marching procedure 
and a finite volume approach. The numerical domain is discretized using C grids, for a good 
description of the leading edge geometry (fig l)- The unknowns are located at the vertices of 
the mesh cells, in such a way to avoid extrapolation towards the blade wall. The code can 
handle open trailing edges, in order to allow additional degrees of freedom in the geometrical 
definition of the blades. This makes the problem of solution existence easier to solve and 
allows a sufficient blade thickness to contain the boundary layer. 

The equations are integrated in time using a Runge-Kutta first order accurate scheme, 
with local time-stepping, enthalpy damping and implicit residual averaging to accelerate the 
convergence. A detailed description of the solver may be found in Leonard (1990). 

CALCULATION OF THE UNKNOWNS ON THE BLADE WALL 

The method developed by the authors is an ’’iterative inverse method”, in which the filial 
geometry is the result of the flow calculation, imposing the required Mach number distribution 
on the blade wall. It has to be iterative since the location of this boundary is part of the 
solution, approximated at the beginning of the design procedure by any convenient initial 
geometry. There is no reason that the flow remains tangent to this geometry during the 
calculation, except in two cases, when the blade ”is” the solution of the problem or when the 
blade wall is modified in order to respect the slip condition, as the time marching procedure 
iterates to the steady state. 

Methods based on the second case have been proposed by Meauze (1982), Giles and Drela 
(1987), and Zannetti et al. (1984). This approach has not been considered here since a 
minimum of successive blade modifications and corresponding mesh generations is desired. 
The blade wall must therefore be treated as permeable to the flow field. After convergence of 
the time marching procedure, the flow calculation results in a distribution of a normal velocity 
component on the blade wall that is used to modify the geometry. 

The calculation of the unknowns at a boundary is dominated by the mathematical nature 
and the physical properties of the system of equations. As the Euler unsteady equations 
are hyperbolic, the solution can be constructed, at any location in the calculation domain 
(including the boundaries) using the information propagating in directions perpendicular to 
characteristic surfaces. The eigenvalues_of tlm Jacobian matrices of the Euler system, projected 
in a considered direction n, are Vn, Vn, Vn a. and Vn. a, and define the propagation 
speeds in that direction. If the vector n is chosen perpendicular and entering the blade wall, a 
positive speed means that the information is propagated on the wave front, in the n direction, 
from the inside of the calculation domain to the outside, and is therefore available to calculate 
the value of the unknowns at this point of the blade wall. On the other hand, a negative speed 
means that the information comes from the outside of the numerical domain and propagates 
towards the inside. This entering information has to be provided by a boundary condition at 
I the boundary point. 
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If the slip condition is imposed on the blade wall ( Vn = 0) only the speed Vn-a is negative 
and therefore only one boundary condition must be imposed, i.e. the velocity direction at that 
point. This shows that the slip condition can not be imposed together with the Mach number 
value, at least for a fixed blade wall. 

On the other hand, if the static pressure p (or the Mach number) is imposed on the blade, 
a velocity component normal to the blade can appear and, depending on its sign, 1 or 3 
eigenvalues will be negative and 0 or 2 additional conditions must be imposed. The sign of 
this normal velocity component can be determined as a function of the imposed static pressure, 
using the compatibility relation corresponding to the only eigenvalue which is always positive 
On + a). 

If the normal velocity is positive, one boundary condition (the required static pressure) 
must be imposed, since only one eigenvalue ( V n - a) is negative. The additional informa- 
tion necessary to calculate all the unknowns at the boundary can be provided by the two 
compatibility relations corresponding to V n and V'„ since they are positive. 

If the normal velocity is negative, two additional boundary conditions must be imposed. 
The best results have been obtained by imposing the total pressure and total temperature at 
that point. Imposing the latter does not give any special problem, since in a blade-to-blade 
calculation it is supposed to remain equal to the total temperature at the inlet. Imposing 
the total pressure is not so straightforward because of numerical dissipation. This problem is 
solved by imposing the value of the total pressure from the previous time level in such a way 
that the total pressure can adapt to the new flow field. This is important when a shock-free 
design is performed starting from a blade for which a shock was present in the original flow 
field, since in this case the initial and final total pressure distributions on the blade wall can 
be very different from each other. A detailed derivation of the compatibility relations can be 
found in Leonard (1990). 

MODIFICATION OF THE GEOMETRY 

A new geometry must be found since the initial shape no longer corresponds to a streamline. 
The modification algorithm is based on a transpiration model and calculates the position of 
the new streamlines using the velocity component normal to the initial blade (Leonard, 1990). 
The modification starts at the stagnation point, and is performed separately for the pressure 
side and the suction side. The new suction and pressure sides are defined as streamlines of the 
flow satisfying the Euler equations, and therefore can not cross each other. This guarantees a 
blade with positive thickness if the numerical integration procedure and the normal velocity 
calculation are sufficiently accurate. 

RESULTS 

The first example illustrates the accuracy of the method for shock-free transonic flows by 
applying it to the supercritical compressor blade designed by Sanz (1984) with an odograph 
method, and proposed as a test case for inviscid calculation methods in AGARD-AR-275 (fig. 
2a). Analysis of the flow with the present method shows discrepancies on the suction side 
Mach number distribution (fig. 2b) similar to the ones observed by Denton (1983). 

The geometry calculated by Sanz has been redesigned using the present method in order 
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to obtain the shock-free Mach number distribution imposed by Sanz as the input data of his 
design (fig. 2c). Only one modification of the geometry has been necessary to obtain good 
agreement (fig. 2d). The difference between the initial geometry designed by Sanz and the 
one designed with the present method is very small. This example suggests that the original 
geometry defined by Sanz may not be shock-free, and illustrates how supersonic flows are very 
sensitive to geometry changes. 

A second example illustrates the design of a shock-free compressor blade, using a NACA-65 
(12A 2 Igb) 1 0 as an initial geometry. This blade is not suited to transonic flows, and a strong 
shock is present in the flow field. Therefore large geometry modifications are expected. The 
flow conditions are: M\ ~ 0.8, p® =1.33 bar, X\° = 341.5 K, (3\ = 45 deg, A/ 2 = 0.5. The 
cascade geometry is defined by a stagger angle of 31 deg and a solidity of 1. 

In a first design, only the suction side Mach number distribution has been modified. The 
initial distribution is compared to the shock-free required distribution in figure 3a. Good 
agreement between the calculated and the imposed Mach number distributions is obtained 
after 4 blade modifications (fig. 3b). The final blade is compared to the NACA-65 blade in 
figure 3c. One observes a thick leading edge, due to the velocity peak in the pressure side 
leading edge region. This is not desirable because it leads to strong diffusion and subsequent 
flow separation along the pressure side, as predicted by a boundary layer calculation. 

A second design has been performed, starting also from the NACA-65 blade, but by modi- 
fying both the pressure and suction side Mach number distributions (fig. 4a). Decreasing the 
pressure side velocity in the leading edge region results in a lower average velocity, and in a 
smaller leading edge thickness because continuity requires a smaller blade blockage. Conver- 
gence to the required distribution is obtained after 3 modifications of the geometry (fig. 4b). 
The initial and final geometries are compared in figure 4c. One can observe a thinner leading 
edge and a shift of the maximum thickness location towards the middle of the blade. 

Blade shapes designed by inviscid methods include the boundary layer blockage on the 
pressure and suction sides. The physical blade geometry can be obtained by subtracting the 
boundary layer displacement thickness from the so-called ’’inviscid” geometry. The minimum 
thickness of the ’’inviscid” blade, required to contain the boundary layer, can be calculated as 
a function of the target velocity distribution before the design procedure is started. 

The analysis of the boundary layer for the prescribed Mach number distribution shown on 
figure 4b indicates that the boundary layer thickness at the trailing edge is of the order of 5% 
of the chord length, which is larger than the total trailing edge thickness of the blade shown 
on figure 4c and makes this blade unphysical. Increasing the trailing edge thickness is possible 
by increasing both the suction and pressure side Mach number distributions in the trailing 
edge region by the same amount (fig. 5a). The circulation is unchanged, resulting in the same 
turning of the flow, but the bade thickness must increase to satisfy continuity. The redesigned 
geometry is compared to the previous design in figure 5b and shows a larger trailing edge 
thickness capable of enclosing the boundary layer and the mechanical thickness. 

An off-design analysis of the second blade has been performed with the same solver, chang- 
ing the incidence by i 2 degrees (fig. 6a and 6b). One can observe that the shock reappears. 
Although the flow field is no longer shock-free, the off-design behaviour of the new blade is 
better than that of the initial geometry. 
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! The third example illustrates the redesign of a transonic turbine blade. The starting geom- 
etry is taken from the workshop VKI-LS 82-05 (Arts, 1982). The flow conditions are: p° = 1 
bar, Tj° = 278 K, /?i = 0 deg, M 2 = 1.1. The cascade geometry is defined by a stagger angle 
of -60 deg and a solidity of 1.25. 

The imposed shock-free Mach number distribution assures a monotonically increasing ve- 
locity on the suction side (fig. 7a). Two modifications of the blade geometry are sufficient to 
give good agreement between the calculated and the required Mach number distributions (fig. 
7b). The original and final geometries are compared in figure (7c). Off-design distributions 
are shown in figures 7d for an exit Mach number of 1.05 and 1.15 instead of 1.1. 

The number of grid nodes used in the above examples ranges from 161 x 15 to 199 x 15. The 
typical amount of CPU time for one blade modification is 15 minutes on an ALLIANT FX/8 
computer with 5 processors. 

CONCLUSIONS 

The present method has been successfully used to design shock-free transonic blades. It 
provides in few iterations results that could not be achieved using traditional direct methods 
and empirical blade modifications. 

The method combines the advantages of a pure inverse method, since the Mach number 
distribution can be imposed on the blade wall, and the advantages of a direct method, allowing 
good control of the geometrical parameters. 

It has been shown how modifications of the required Mach number distribution influence 
the blade geometry. Special attention was given to design trailing edges of sufficient thickness 
to enclose the boundary layer blockage. 

Off design analysis of designed geometries illustrate the difficulty of optimizing for more 
than one operating point. 
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Fig. 1 C-grid discretization 
for a turbine blade 
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Fig. 4a Initial ( + ) and required (Q) 
Mach number distributions 



Fig. 4b Calculated ( + ) and required (O) 
Mach number distributions 
after 3 modifications 



Fig. 4c Original and final 
geometries 
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Fig. 5a Initial ( + ) and required (Q) 

Mach number distributions for 
increased blade thickness 
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Fig. 6a Off design 

(-2 deg incidence) 


Fig. 6b Off design 

( + 2 deg incidence) 
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Fig. 7a Initial ( + ) and required (Q) 
Mach number distributions 


Fig. 7b Calculated ( -i-) and required (0) 
Mach number distributions 
after 2 modifications 




Fig. 7c Original and final 
geometries 


Fig. 7d Off design for M 2 = 105 ( + ) 
and M 2 = 1.15 (O) 









